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Abstract: The free energy of magnetization is expressed as a fourth degree 
polynomial in the magnetization with a linear term due to an 
external field and ctuadratic and q.uartic terms due to Fermi j 
exchange, magnet ostatic, etCoo energies as well as to entropy- 
effect So 

Below the Gurie temperature the free energy function 
has three real-4xtrema, namely two minima and a maximum, if the 
absolute magnitude of the external field lies below a critical 
value* If the absolute magnitude of the external field exceeds 
this critical value, there is only one extremum, namely a mini- 
mumo When the absolute value of the field equals the critical 
value, there is a minimum and a point of inflectiono 

It is assumed that the magnetization of a system is 
always such as to minimize the free energy, but that the system 
is not always capable of settling in the lowest relative minimtm 
of the free energy. On the basis of this model, the critical 
field and the correspoiilding magnetization are calculated as 
functions of parameters peculiar to the material and temperature, 
and experimental data obtained from hysteresis loops at various 
temperatures are presentedo 

Calculation of magnetization as a function of the 
external field results in a universal expression for the hysteresis 
loopj vrhich is compared with some experimentally observed loops. 

Introduction 

It can be shown that the internal energy of magnetization in the 
absence of an external field can be expanded in an even series in the magnetization 

Xpt (A ^ B Hh 0) M^ ♦ in ♦ fto». 



St 



la Slater; •quantum Theory of Matter *• Mc&raw-£ill Book Co«« Inc. 1 Edo (1951) 
Appendix 22^ p*» 5154, 
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^2 V 

where: AM is a term due to the Permi energy; A> 

BIT is a term due to the exchange energy; B > 

2 

CM is a term due to the magnet ostrictive, magnetostatic , etc, energies 

k 
DM includes higher order terms in all these energies. 

Similarly, the entropy can be expanded as follows: 

S = -PM^ - QM^ - ..... P >0, <%>0 
The free energy is then given by: 

P = B - TS = (A - B ♦ C 4- PT)M^ Hh (B 4- qT)M^ (Equ. 1) 

which has extrema for the values 



«°o.t\/ -"/(g;§ ^ 



PT 



A ■> B 4> C 
When 3? ^ • — ^^-5 — — there is only one realminimum and the material is para- 
magnetic. 

A •> B 4" C 
When T ^ — ^ • — there are two real minima separated by a maximum 

at M = 0» The system is in this case ecLually likely to assume a positive as 
a negative magnetization. In fact some regions may become positively magnetized, 
others negatively, so that the total magnetization is zero. Seversing the 
magnetization of a region requires the energy needed to pass over the maximuta 
at M = 0. Since such reversal is improbable » the magnetization of anyone region 
remains constant during an observation, so that domain patterns are obtained 
for ferromagnetic, but not for paramagnetic materials. Por the latter, there 
is no energy barrier; therefore, anyone region is equally likely to have a 
positive as a negative magnetization with reversal of magnetization occurring 
all the time. Daring the period of an observation anyone region in a paramag- 
netic material has therefore an average magnetization zero. 

The transition from ferro- to paramagnetism occurs at a temperature 
T = T given by: 

^c P 



This is by definition the Curie temperature. Substitution into 

equation 1 produces: 

P = -P (T^ - T)M^ + (D + QT)M^ (Equ. 2) 

c 
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ifhen an external f ield H is applied to the system, the free energy- 
is decreased by an amount HM« so that 

2 4 

F = -HM - P(T^ - 5?)M + (D + WU CBq.u. 3) 

S'igu.res la, lb, Ic, Id, le. If, and Ig show the free energy plotted vs. mag- 
netization for various values of the field strength. 

In previously unmagaeti zed mp.terials, the total volume of the domains 
with positive magnetization equals that of the domains with negative magnetiza- 
tion, the magnitude of each magnetization being given by the position of the 
minima in Figure Id, When a positive field is imposed, the free energy changed 
as indicated by Figures Ic, lb, and la. The positions of the minima are changed 
by the application of a field; a positive field causes both minima to shift 
in the positive direction,, so that the volume of positively magnetized domains 
now exceeds that of negatively magnetized domains. Some motion of domain 
walls is thus required to maintain the system at minimum, free energy. Figure 
lb shows the free energy c\arve for the critical free strength H » H where 
there is only one minimum as well as a point of inflection. When H^ H , there 
is no energy barrier keeping the system from settling in the lower (now only) 
relative minimum. Therefore when the field strength increases from zero, 
a cataclysmic reversal of the domains with magnetization opposite in direction 
to the external field occurs when the field strength equals H - Further 
increase in field strength simply shifts the position of the minimum to the 
right, so that the average magnetization increases toward saturation, 

A subsequent decrease in field strength decreases the average 

magnetization, 1 When H = H is reached, the free-energy curve again has the 

c 

form shown in Figure lb, amd further decrease in field strength produces 
successively the curves in Figures Ic, Id, and le. The magnetization of 
the system is now given by the position of the right-hand minimtua, which is 
the lower of the two relative minima, when H> 0, and the higher when H<^0, 
When E < there is, therefore, a tendency for the system to reverse its 
magnetization, which is prevented by the energy barrier. As shown ia Figure If, 
this barrier "vanishes when H = -B , so that another cataclysmic reversal then 
takes place. ILlternating the magnetic field with an amplitude greater than 
[hi therefore produces a nonlinear variation in magnetization such as that 
described by the hysteresis loop. The upper branch of the loop corresponds 
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to the right-hand minimum, hence to a decreasing field; the lower hranch 
corresponds to the left-hand minimum and increasing field. When |H|>Ih I, 
the magnetization curve is single "valued because the free energy has only 
a single minimum. 

Computation of the Critical Field and the Corresponding Magnetization 
According to what has preceeded> the hysteresis curve can be 
calculated on the basis of the assumption that the magnetization has a value 
that minimizes the free energy. The vialues of magnetization that minimize 
the free energy are obtained by differentiating the free energy as given in 
equation 3 with respect to the magnetization and equating the result to zero: 

2 
The discriminant Of this equation is: 



A- fHTr. - T) l3 27 r H - ] 
2 L B + QT J "1^ Lb + QT J 



A double root, i.e. a point of inflection in the free energy curve, as shown 

in Pignres lb and If, occurs when this discriminant equals zero. Therefore 

the critical value of the field strength, H , is found by equating the discrim- 

c 

inant to zero. 

.. H , |^/]Cp(??r 7 ^)1 (Equ 5) 

• • ^0 /B-»- QTVa V-BJqu-. ^i 

Substituting this value for H in equation 4 and solving the latter 
equation for M produces the value of the ma^etization just after the cata- 
clysmic reversal has occurred, which is called M : 

% DTqf "13 mn^ — J ^^"^^ ^^ 

ffrom equations 5 aJid. 6 can be derived the following two relations: 

IT = (2/3] P(T, - f ) («iu. 7) 



'"^ la. = D + qj (»iu. 8) 



«o 



3 
By plotting experimentally obtained values of Hq/Mq and Hq/Mq'^ vs. T, one 

can calculate the coefficients P, B, and Q (See "Results and Comparison with 

Experimental Data). 



2, «JHevr Pirst Course in thi Theory of; Equations* L. 1. Bickson, J. Wiley & 
Sons, p. 46, 1939. 
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Computatioa of the Masaetization Curve 

To find the magaetization for any value of the field strength, 

it is necessary to solve equation 4 for all values of H« Equation 4 can he 

rewritten by substituting for P(T -T) and (D + QT) the values given by equa- 

c 

tions 7 and 8: 

,3 3 o 2 „ 1 Ur^ 



M-^ -^ M 



M - -i--ifi- H 



c " - -r ^c 
^1) -4 <r> -T'4> = ° ^^^- 9) 

c ' c c 
The variables M/M and H/H express the magnetization and critical field 

strength as dimensionless quantities in units of the critical values M and E . 

This is just ifhat is done experimentally when the oscilloscope picture is 

normalized, Equation 9 can therefore be said to represent a normalized 

hysteresis loop; the normalizing coefficients H and M depend on the tern- 

perature and the nature of the medi\a^ as expressed by equations 5t 6, ?, and 8, 

but equation 9 is universal, i.e. independent of the medium and of the 

temperature. 

Equation 9 must be solved separately for the case | H(£|H J and for 

the case |H|>|H^|. For the former case, set H/H = cos 3A. Then the three 

real solutions are: 

^ = cos A, cos (A + 120°), cos (A ■*- 240°) , (Bqu. 10) 

where the first two solutions correspond to free energy minima, the latter 
to the energy maximum. Since equation 9 contains no quadratic term, the siam 
of its three solutions equals zero, so that the energy maximum occurs at a 
value of the magaetization equal to the difference in absolute values of the 
magaetization corresponding to the energy minima. 

When |h|>|hJ , M/M^ is given by equation 11. 



I- = ^ >1 f- \l Ii-Sl l^{^^ 4.7 1. \/l-(|fi)^' (Equ. 11) 
c c L "^ 

and hence approaches a cubic parabola asymptotically, , ifhen H = i H , 



M » ± M according to both eqiaations 10 and 11, 



♦ This can be proven by direct substitution into equation 9, remembering that 
cos"^ X =/l/4j(cos 3x4-3 cos x). 
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Res-ults and Comparisons with Experiments 

Figure 2 shows a plot of H^M vs. (3? -T) using experimental values 

reported "by Channing Morrison for a Mn-Kg ferrite (MJ-lllS, I'-.259). Hear 

the Curie point the horizontal and vertical gstin are very large, so that 

failure to reach saturation is overstressed in this region. This results 

in too large values of H /M and accoimts for the tendency of experimental 

c c 

points to veer away from the origin. 

Jo ,■ 3 

Prom Figure 2 it is found that P = 8.3^ x 10 erg /cm-^/ degree. 

gauss'^ 

It is hoped that during the .summer of 1953 sufficiently accurate measurements 

3 
will he made so that H /M can also he plotted in order to botain I) and C^. 

c c 

It will he interesting to compare results of temperature dependent measurements 
for various materials. 

In Figure 3 is drawn the hysteresis loop obtained by solving 
eq^uation 9 , and for comparison Figure ^ shows a loop obtained experimentally 
with molybdemam permalloy. 
Conclusions 

It has been shown that simple fundamental models for magnetism 
lead directly to a hysteresis loop which resembles many experimentally 
observed hysteresis loops. The calculated loop has the characteristic 
■Sloping shoulders,* which have heretofore been explained by a rotation of 
domains. The preseni model neglects such rotation, and yet the sloping shoulder 
is apparent. 

The very simple model employed has led to the universal equation 9 
in much the same way as a simple kinetic gas theory led to the universal 
Boyle' s law. Deviations from the universal law depend on the individual 
materials concerned, Some of the causes of deviations from the universal 
hysteresis loop will now be discussed. In the first place, the experimentally 
observed loop represents dft/dt vs. I, and the finite size of the sample causes 
the sides of the loop to slope. Secondly, it has been assumed that reversal 
of magnetization occurs only when the energy barrier vanishes. Actually a 
■tunnel effect* allows reversal through the barrier, the speed (hence amount) 
of reversal being proportional to the negative magnitude of the barrier 
expressed in units of kT. Thus the sides of the loop do slope somewhat be- 
cause the reversal can start somewhat below the critical field strength. 
Minor loops are obtained when the amplitude of the alternating field strength 



3. Ohanning Morrison, :|LiJ91, "Hysteresis Loop Characteristics of MF-1118 
for Different Temperatures," Digital Computer laboratory, MIT, 
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is slightly less than the critical value; partial reversal then takes place 
through the harrier, pat a good deal of the material remains captured hehind 
the Ijarrier, so that the resulting magaetization is less than on the satura- 
tion loop, 

Figure 6 shows a "DC" loop. Ifhen the system is reversing its 
magnetization (|H| £ JH |) the magnitude field strength is suddenly decreased 
in magnitude p so that an energy harrier is reintroduced^ Ihe material that 
had not yet reversed its magnetization is then caught again, Imt the material 
that was already reversed is now in the lower relative minimum, and hence will 
not reverse again. Therefore the magnetization remains essentially constant 
\mtil the magnitude of the field strength is once more increased to the critical 
value, when more reversal can take place, 

Thirdly, only two directions of magnetization have "been assumed, 
namely positive and negative. Actually a range of 360 with respect to the 
direction of the external field is possible, and Figure 1 should he three 
dimensional with a free energy surface analogous to lyring surfaces in reaction 
kinetics. The system might reverse its magnetization by ^oing around the 
energy barrier rather than over the topj this would amount to rotation of the 
magnetic vector. lor directions of easy magnetizations the free energy surface 
would have a pit surrounded by barriers; when the material is isotropic, 
Figure Id would only need to be rotated about its axis of symmetry to produce 
the free energy surface. Thus the hysteresis loop could be determined by 
the motion of a particle in a complicated system of mountains, crevices and 
passes, with deviations from a straight path representing rotation. The 
possibility of rotation would make the loop less sa^jar© than the one shown 
in Figure 3« Further it has been assumed that the only force opposing re- 
versal of magnetization is that due to the free energy barrier, and that the 
reversal occurs instantaneously. Frictional forces, i,e, , forces that are 
not simply derived from a potential as assumed in eci_uation 1, have not been 
taken into accoun|;. 

Finally, it is assumed that the system can always be in a state 
of lowest free energy. Actually * super saturation* may occur, for the 
mechansim of reversal of magnetization requires formation of domains of 
reverse magnetization just as precipitation form solution requires nuclei 
of condensation to keep the solution from being supersaturated. The 



^. J. B, G-oodenough, E-.532, •Nucleation of l>omains of Reverse Magnetization 
& Switching Characteristics of Magnetic Materials,* Digital Computer 
Laboratory, MIT, 
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third aad fourth factors tend to make the loop more sq,\iare than that shown 
in Pigare 3* Figure 5 shows an experimentally observed loop with »' double 
shoulders,^^ The shoulder can be explained as follows: The observed hystere- 
sis loop is a combination of two loops. The outer loop represents one in which 
reversal of magnetization is inhibited by the lack of domains of reverse 
magnetization. The inner loop represents the free energy or equilibrium 
hysteresis loop of Figure 3, . When the absolute value of the field reaches a 
critical value , domains of reverse magnetization ere formed, so that the 
magnetization drops down to the equilibrium loop. 



ALL/jk 
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FIG. 1 

FREE ENERGY OF MAGNETIZATION IN AN EXTERNAL FIELD 
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REDUCED HYSTERESIS LOOP 
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FIG. 5 

PORTION OF A DC HYSTERESIS LOOP WITH 

MULTIPLE SHOULDER 
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FIG. 6 

HYSTERESIS LOOP OF PERMINVAR 

(REPRINTED FROM J. APPL. PHYS. 23, 318 
WITH PERMISSION OF THE AUTHORS, M.J.WILLIAMS AND M. GOERTZ) 
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